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Abstract 

A computational heat transfer methodology was developed to study the dual-engine linear 
aerospike plume induced base-heating environment during one power-pack out, in ascent flight. 
One power-pack out results in reduction of power levels for both engines. That, in turn, reduces the 
amount of base-bleed and changes the distribution of base-bleed on the two pillows. Hence, the 
concern of increased base-heating during power-pack out. The thermo-flowfield of the entire 
vehicle was computed. The computational methodology for the convective heating is based on a 
three-dimensional, finite-volume, viscous, chemically reacting, and pressure-based computational 
fluid dynamics formulation. The computational methodology for the radiative heating is based on a 
three-dimensional, finite-volume, and spectral-line-based weighted-sum-of-gray-gases absorption 
computational radiation heat transfer formulation. A separate radiation model was used for 
diagnostic purposes. The computational methodology was systematically benchmarked. In this 
study, near-base radiative heat fluxes were computed and they compared well with those measured 
from an installed linear aerospike engine tests. The base-heating environment of 18 trajectory 
points selected from three power-pack out ascent scenarios was computed and is presented here. 
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The power-pack out condition has the most impact on convective base-heating when it happens 
early in flight. The source of its impact comes from the asymmetric and reduced base-bleed. 

Nomenclature 

C,AAA= turbulence modeling constants, 1.15, 1.9, 0.25, and 0.09. 

G = geometrical matrices 

H = total enthalpy 

h = static enthalpy or altitude, km 

I = radiative intensity 

J = Jacobian of coordinate transformation 

K = forward rate constant 

k = turbulent kinetic energy 

M = Mach number 

N = total number of chemical species 

P = pressure 

Pr = Prandtl number 

Q = heat flux, kW/m 2 

q =l,u,v,w,H,k,e,orpi 

R = recovery factor 

r = location coordinate 

S 4 = source term for equation q 

T* = law-of-the-wall temperature 


t 


= time, s 


u 


= volume-weighted contravariant velocity 
= mean velocities in three directions 


u, v, w 

u t = wall friction velocity 

u + = law-of-the-wall velocity, (u/u x ) 

x,y,z = coordinate or distance 

y + = law-of-the-wall distance, (y p u t p/|i) 

e = turbulent kinetic energy dissipation rate or wall emissivity 

k = absorption coefficient 

H = effective viscosity, (pi + m) 

% = computational coordinates 

IT = turbulent kinetic energy production 

p = density 

<j q = turbulence modeling constants 

<p = energy dissipation function 

Q = direction vector 

co = chemical species production rate 

Subscripts 

a = ambient 

b = black body or base 

c = convective or center 

EA = Engine A or “on” engine 


EB = Engine B or “off” engine 

1 = laminar flow 

p = off-wall (wall function) point 

r = radiative 

t = turbulent flow 

w = wall surface 

«= = freestream 

0 = reference 


Introduction 

X-33, a half-scale prototype of a proposed commercial, reusable launch vehicle (RLV), was 
designed to demonstrate new, reusable single-stage-to-orbit technologies. It was fueled by liquid 
hydrogen and liquid oxygen and featured a lifting body configuration coupled with two 
integrated, linear aerospike rocket engines. One goal of the project was to provide safe, reliable 
and affordable access to space. 

It is well known that aerospike engines have the potential advantage of adjusting themselves 
to perform with maximum efficiency at all altitudes. 1,2 However, the aerospike plume- induced, 
base-heating environments must be understood fully in order to satisfy safety and reliability 
requirements. One scenario of particular interest is the base-heating environment during a 
potential power-pack out (PPO) - the loss of one of the two integrated turbopump sets - resulting 
in reduced power levels for both engines. This in turn reduces the amount of base-bleed which 
protects the plug-base components such as the pillows and flex seals (see Fig. 1). The reduction of 



base-bleed during PPO is a concern in light of observed heat-induced flex seal damage, seen after 
a single-engine, hot-fire test at reduced power levels. Furthermore, the plumbing is such that for 
PPO operation the “on” engine-side of the pillow receives about 75% of the remaining base- 
bleed, whereas only 25% goes to the “off’ engine side. This creates a concern for the integrity of 
the plug-base components and the need to understand the physics and the impact of PPO on the 
ascent base-heating environment. Note that subscale testing is not useful since combustion 
effects such as radiative heating are not scalable. In addition, full scale testing is not feasible 
(since the risk of losing an engine is not acceptable) and cannot simulate altitude effect. Using 
anchored numerical tools to computationally simulate the base-heating environment during PPO 
is therefore a reasonable strategy to address these concern and needs. 

In the following sections, the computational methodology and the benchmarks performed are 
delineated first. The PPO trajectory and the run matrix are then discussed. And finally, the 
computational results showing the effect of base-bleed on base-heating environment during PPO 
are presented. 


Solution Methodology 
Computational Grid Generation 

Previous studies 3,4 have shown that the X-33 base-flowfield is three-dimensional due to the 
linear aerospike layout and the lifting body configuration. And, only half of the domain was 
computed since flow symmetry can be assumed in normal operating conditions. In those studies, 
the computational grid was gradually improved through a sequential grid-refinement technique to 
an eventual density of 2,217,444 points. 3 Under PPO situations the flow is asymmetric, hence 
the flow over the entire vehicle is solved. The grid density is hence doubled to 4,434,888 points 



by mirror-imaging the original half-vehicle grid using software package GRIDGEN. 5 Figure 1 
shows a layout of the surface computational grid with component definitions. Notice the plug- 
base flex-seal region where severe base-heating is expected after PPO. An identical grid is used 
in the radiation calculation. The total grid-density is increased to 4,724,592 points however, due 
to the accounting system in the radiation formulation. 

Flow and Radiation Computation 

The effect of base-bleed on the X-33 base-heating environment is investigated with two 
computational tools: the Finite-Difference Navier-Stokes (FDNS) CFD code 6 for the convective 
heating and the General Radiation Solution Program (GRASP) 7 for the radiative heating. These 
tools were developed at Marshall Space Flight Center (MSFC) and are continuously being 
improved by MSFC personnel and its supporting contractors. Systematic and rigorous 
benchmark studies have been performed for base flow and heat transfer applications. For 
example, FDNS has been validated for convective heat transfer inside rocket thrust chambers 8 
and coolant channels 9 , for base-pressure characterization of a four-engine clustered nozzle 
configuration, 1011 for Delta Clipper-Experimental (DC-X) base-drag induced by the engine 
exhaust during cold-flow and flight tests 12 , and for DC-X convective base-heat flux during 
landing 13 . GRASP has been benchmarked for DC-X radiative base-heat flux during landing 13 . 
In this study, FDNS and GRASP calculations are conducted sequentially in order to save 
computational memory. The CFD computed pressure, H 2 0 concentration, and temperature fields 
are inputs for the radiation calculations. The solution algorithm for the flow and radiation 
computation is summarized in the following. 





Flow Solution and Convective Heat Transfer 


FDNS solves a general curvilinear coordinate, chemically reacting, viscous thermo-flowfield 
with Reynolds-averaged transport equations. A generalized form of these equations is given by 


dpq j\rPUq^lO q )3{dqldCl\ S q 
Jdt d£ + j 


A pressure-based predictor-plus-multicorrector solution method is applied. 14 The basic idea is to 
perform correction for the pressure and velocity fields by solving for a pressure correction so that 
the velocity-pressure coupling is enforced, based on the continuity constraint. A second-order 
central-difference scheme is employed to discretize the diffusion fluxes and source terms of the 
governing equations. For the convective terms, a second-order total-variation-diminishing 
difference scheme is used in this effort. 

An extended k-e turbulence model 15 is used to describe the turbulence, pt = pC 4 k 2 /e is 
defined as the turbulence eddy viscosity. Turbulence modeling constants a q and source terms S q 
of the transport equations are given in Table 1. These turbulence modeling constants have been 
used extensively for combustion-driven and base flows 8 ' 14 ’ 16 , while o k and a e are taken from the 
turbulence closure. 15 A 7-species, 9-reaction detailed mechanism 16 (Table 2) is used to describe 
the finite-rate hydrogen-oxygen afterburning chemical kinetics. The seven species are H 2 , O 2 , 
H 2 O, O, H, OH, and N 2 , where H 2 O is the major radiating medium. Although it has been shown 


that th e- 


power flight, 3 the finite-rate chemistry is turned on for all the trajectory points in this study due to 
considerations of reduced power level and base-bleed during PPO. 



Table 1 o q and S q of the transport equations 


q 

CT q 

Sq 

i 

1.00 

0 

U 

1.00 

-Px+V[p(Uj) x ]-(2/3)(pVuj) x 

V 

1.00 

-Py+V[p(Uj) y ]-(2/3)(|xVUj) y 

w 

1.00 

-Pz+V[jx(Uj) z ]-(2/3)(p,Vuj) 2 

H 

0.90 

DP/DT-np 

k 

0.89 

p(n-e) 

e 

1.15 

p(e/k) { [Ci+C 3 (TI/e)3n-C 2 e } 

P« 

0.90 

(Bi, i = 1,...,N 


Table 2 Hydrogen/Oxygen combustion kinetics mechanism. Kf = AT® 


Reaction 

A 

B 

E/R 

h 2 +o 2 =oh+oh 

1.700E13 

0 

2.407E4 

oh+h 2 =h 2 o+h 

2.190E13 

0 

2.590E3 

oh+oh=o+h 2 o 

6.023E12 

0 

5.500E2 

o+h 2 =h+oh 

1.800E10 

1.0 

4.480E3 

h+o 2 =o+oh 

1.220E17 

-0.91 

8.369E3 


M + 0 + 0 = 0 2 + M 2.550E18 -1.0 5.939E4 


-E/RT 


M + H + H = H 2 + M 5.000E15 0 


0 



M + H + 0H = H 2 0 + M 8.400E21 -2.0 


0 


* M stands for third-body collision partner 

A modified wall function approach is employed to provide wall boundary layer solutions that 
are less sensitive to the near-wall grid spacing. Consequently, the model has combined the 
advantages of both the integrated-to-the-wall approach and the conventional law-of-the-wall 
approach by incorporating a complete velocity profile 17 given by 

. \ 4.02 ( .2 . \0.T9 

y + +ll) /[ / -737/ +833 1 
+5.63tan‘ , (<>.12/ -0.44l)- 331 

and a universal temperature profile 18 given by 

T + = U + + 1230*^ 068 - 1) (3) 

The convective heat transfer from a hot boundary layer to a cooler wall follows the modified 
Newtonian law 13 

Q cw = (pu/T)[h w - h p - R(u p 2 /2)] (4) 

where R = Pq 1/2 if y + < 1 1.63 and R = Prj 1/3 if y + > 1 1.63, and y + = 1 1.63 is the thickness of the 
viscous sublayer. A constant Pq of air is used in this study, since parametric studies performed in 






Ref. 13 show that the H 2 /O 2 plume induced base heat flux is not sensitive to a multicomponent 
variable Pri. 

Radiative Heat Transfer 

GRASP analyzes the radiative field by solving the general curvilinear coordinate radiative 
transfer equation with a finite-volume method (FVM) formulation: 19 

(QV)/(r,Q) = -K/(r,£2)+K/ 4 (r) (5) 

The term on the left-hand side represents the gradient of the intensity in the direction of Q. The 
two terms on the right-hand side represent the changes in intensity due to absorption and 
emission. The wall boundary is assumed gray while emitting and reflecting diffusely, and the 
radiative wall boundary condition is given by 

/(r„,Q+) = eI b (r w ) + £zfl f /(r M> ,Q-)|i,.Q-| < «- (6) 

It J 
n.Q~<0 

with 

= J/(r*,Q")|n.£r|<jQ - (7) 

n.Q~<0 

where £2 + and £2' denote the leaving and arriving radiative intensity directions, respectively. The 



20-band spectral-line weighted-sum-of-gray-gases model 7 is used to calculate the total emissivity 
and absorptivity of the radiating medium. Following the ray-dependency test performed in Ref. 
13, the FVM 6x4 option - six control angles in the polar direction and four in the azimuthal 
direction - is deemed adequate and is used in this effort. 

Boundary and Initial Conditions 

Hie outer boundaries of the computational domain have the following conditions: the fixed 
total-condition (free-stream) and the flow exit plane. A no-slip wall is specified for the body 
surface. A fixed (ambient) static pressure is imposed on the exit plane and on a point far away 
from the action area (one grid point off the freestream boundary), in order to obtain a unique 
solution for the desired altitude. The fixed inlet boundary condition is applied to the exit plane of 
the 40 thrusters where the flow properties were obtained from a chemical equilibrium analysis 
using the CEC code. 20 

Base-bleed is accomplished in the actual engines by flowing a mixture of hydrogen gas and 
steam from thousands of orifices drilled on the pillows. Because it is computationally 
prohibitive to resolve all the individual orifices in the CFD model, the orifices and the pillow 
surface surrounding them can be lumped together and modeled as one base-bleed inlet boundary. 
In Ref 3 and 4, a fixed mass flow-rate (and energy) boundary condition is used in subsonic 
freestream cases and a fixed inlet boundary condition is applied in supersonic freestream cases. 
This approach however, precludes the calculation of convective heat flux on the pillow since it is 
treated as an inlet and not a wall. In this study, a special base-bleed boundary treatment is 
devised by combining the features of constant mass flow injection normal to the surface and 
those of wall boundary treatment such that the increased convective heating due to reduced base- 



bleed can be estimated for subsonic freestream cases. For supersonic freestream cases, the pillow 
convective heating is negligible and the fixed inlet boundary condition is retained. The basic 
assumptions for the base-bleed boundary treatment are that the injection is always normal to the 
wall surface and the mass flow rate specified is constant throughout the computation. Constant 
injection temperature and species composition are imposed to maintain a fixed inlet enthalpy. 
Since the base-bleed boundary is treated partially as a wall, the wall function approach is applied 
and the wall heat flux can be estimated. Figure 2 shows the computed base-bleed flow vectors at 
40 s into the PPO at launch + 30 s trajectory. Note that the flow vectors (mass flow rate) of the 
“off” engine side are about 1/3 the length of those of the “on’ engine side. 

For convective heat transfer calculations, ambient temperature is prescribed as the forebody 
and aftbody surface temperatures, whereas 540 deg R is specified for all base surfaces per base- 
heating design convention. For radiation calculations, the surface emissivity of the entire vehicle 
is assumed to be 0.7. 13 The engine ramp is actively cooled and the surface temperature 
distribution is prescribed from a separate conjugate heat transfer calculation involving solid walls 
and coolant channel flows. It was found 3 that cowl-base irradiation is more than ten times higher 
if an adiabatic condition is imposed on the ramp surface. Hence, a more accurate boundary 
condition is employed. Notice that the effect of re-irradiation from surface radiation is included 
in all calculations. 


X-33 Specific Computational Benchmarks Performed 
Forebodv and Aftbodv Surface Pressure Coefficients 
The X-33 forebody and aftbody surface pressure coefficients were computed and compared 
reasonably well with those of a 7.75% scaled model cold-flow test. 3 These comparisons indicate 


that the incoming flow environment is adequately simulated for X-33 base-flow development. 

Installed Linear Aerospike Engine Plume Induced Noise 
Hie X-33 CFD plume data base 3 was used in conjunction with an engineering acoustic 
radiation model to predict the acoustic environment 21 of the linear aerospike engine, with that 
measured from a series of hot fire tests. The predicted far-field acoustic spectra, directivity 
pattern, and overall sound pressure levels agreed reasonably well with those of measurement. 
Since the acoustic radiation model uses CFD computed (du/dn) 1 , (k/e), and (kp/<p) for time 
scaling, this benchmark indicates the CFD computed X-33 plume velocity and turbulence fields 
are adequate. 

Base-Pressure Characteristic curves 

The base-pressure characteristic curve is the central base-pressure plotted as a function of 
ambient pressures. The computed base-pressure characteristic curve was compared with that of a 
7.75% scaled model cold flow test 3 and that of a 2.25% scaled model hot flow test. 4 The 
reasonable agreement indicates that the predicted base-pressure characteristic curve is adequate. 

Sea-Level Base Pillow pressures 

The computed sea-level base pillow pressures agreed reasonably well with those of a 2.25% 
scaled model hot-flow test. 4 The predicted and measured 22 pressure levels indicate there is a 


strong aspiration occurring on the base pillow at low altitudes. 




Base Horizontal Centerline Convective Heat Fluxes 


Base horizontal centerline convective heat fluxes were computed and compared with those of a 
2.25% scaled model hot-flow test at sea level, at 23.4 km, and at 23.4 km with fences installed. 4 
The measurements 22 were made on the plug-base, side wall, inner-base, offset outer-base and 
outer-base. Not only is the predicted convective base heating in good agreement with 
experimental data, the linear aerospike base-heating physics such as the formation of the reverse 
jet, lateral wall jet, plume spillage, and aspiration jet are well captured. 

Installed Linear Aerospike Engine Plume Induced Base Radiative Heat Fluxes 

In this study, near base radiative heat fluxes were computed and compared with those 
measured from a series of hot fire tests of a single X-33 XRS-2200 linear aerospike engine at the 
NASA Stennis Space Center in 1999 and 2000. Wide band plume radiation data were collected. 
Plume radiation instrumentation consisted of several 180° hemispherical radiometers. A typical 
view of the engine and plume during test is shown in figure 3. 

Wide-angle (180 deg.) radiometer locations were selected to provide orientations similar 
to that of flight vehicle surfaces and flight test instruments. Most instruments (e.g., radiometer 
R101B, see Fig. 4) were located to the side parallel to the engine long axis in what would be the 
engine base region. One radiometer (R105, see Fig. 5) was in close proximity to the engine near 
the cowl base and thruster area on the ramp side. 

Radiation predictions were made using FDNS flow field solutions with heating rates 
generated by the GASRAD code. The GASRAD code uses a statistical band model for 
exponential line strength distribution with Lorentz/Doppler line shapes for gaseous species. Its 
geometry model provides high fidleity integration over a hemisphere and can accept either 





axisymmetric or three dimensional plumes. GASRAD is particularly well suited for very 
accurately predicting radiation to indiviual geometry points, while GRASP is typically used for 
entire surface predictions due to its relatively high computational efficiency. Details of the 
GASRAD radiation formulations can be found in Ref. 23. 

An FDNS solution was generated at sea level ambient pressure. Ground test radiation 
data are compared to sea level predictions in figures 4 and 5. Predicted levels agreed quite well 
with measured data. Temporal variation of the data is due to engine startup and shutdown 
transients, as well as scheduled power level and fuel mixture ratio changes. 

Some additional insight into modeling accuracy can also be gained from these data. 
Radiometer R101B has a view factor encompassing a large potion of the engine plume. The fact 
that the predicted level agrees well with the data as shown in figure 4 gives us confidence that the 
overall plume flow field is being modeled accurately. Figure 5 shows that data from the cowl 
base region radiometer (R105) is also predicted accurately. This indicates that the local flow 
field just downstream of the engine thrusters is being accurately modeled because the field of 
view of R105 is predominant in this region. The results of these benchmarks validate our 
computational heat transfer model for studying the X-33 aerospike plume induced base-heating 
environment during PPO. 


PPO Trajectories and Run Matrix 

The linear aerospike plume-induced X-33 base-heating environment was computed for trajectory 
points selected from three PPO scenarios: PPO at Launch (L) + 30 s, PPO at L + 60 s, and PPO at L 
+ 120s. The ascent mission was expected to fail if PPO occurs before L + 30 s. Figure 6 shows the 
computational points selected on the thrust, altitude, and Mach number profiles for those three PPO 


trajectories. Six strategic points were selected for each trajectory in order to bracket the base- 
heating environment. Note that one point was chosen before the PPO and another point was chosen 
after the PPO. This is best illustrated in thrust profiles since the biggest power difference occurred 
right after the PPO. To save computational time, several selected points were overlapped on the 
three PPO trajectories. For example, the trajectory points of 0 and 30 s on the PPO at L + 60s 
trajectory are identical to those of 0 and 30 s on the PPO at L + 30 s trajectory. This way, the 
original number of computational points of 18 is reduced to 13. The power levels of the dual 
engines at selected trajectory points were provided. The base-bleed mass flow rates on the two 
pillows were calculated based on the 75/25 split for engine-off cases. Tables 3-5 show the ran 
matrix for the three PPO trajectories. 


Table 3 Run matrix for PPO @ L+30 s trajectory 


Case 

t,s 

M» 

H, ft 

PLea 

PLeb 

O/Fea 

O/Feb 

ntbb.EA 

mbbyEB 

1 

0 

0.00 

0 

0.8187 

0.7971 

5.7225 

5.4298 

9.930 

9.666 

2 

30 

0.44 

9495 

0.9958 

1.0023 

5.6546 

5.6489 

12.088 

12.168 

3 

40 

6.37 

13821 

0.4937 

6.4827 

5.7224 

5.4032 

8.995 

2.919 

■ 

100 

0.21 

28954 

0.4992 

0.4904 

5.5420 

5.6735 

9.119 

2.959 

5 

200 

1.39 

77217 

0.5109 

0.4817 

5.6897 

5.7159 

9.147 

2.968 

6 

280 

5.00 

168498 

0.4571 

0.4043 

6.1580 

5.8515 

7.918 

2.570 


Table 4 Run matrix for PPO @ L+60 s trajectory 


Case 

t,s 

ML 

H, ft 

PLea 

PLeb 

O/Fea 

O/Feb 

Ttlbb % EA 

mbb t EB 



































































1 

0 

0.00 

0 

0.8187 

0.7971 

5.7225 

5.4298 

9.930 

9.666 

2 

30 

0.44 

9495 

0.9958 

1.0023 

5.6546 

5.6489 

12.088 

12.168 

3 

60 

0.74 

28904 

0.7200 

0.7449 

5.6331 

5.7783 

8.7262 

9.0297 

■ 

61 

0.74 

29716 

6.3780 

0.3984 

4.6746 

4.9030 

7.0679 

2.2936 

5 

200 

3.49 

130321 

0.5023 

0.4854 

5.3726 

5.7644 

9.0121 

2.9244 

6 

270 

7.28 

166925 

0.4277 

0.4161 

5.4451 

5.6719 

7.6884 

2.4949 


Table 5 Run matrix for PPO @ Lf 120 s trajectory 


Case 

t,s 

Mco 

H,ft- 

PLea 

PLeb 

O/Fea 

O/Feb 

. 

ntbb.EA 

mbb,EB 

1 

0 

0.00 

0 

0.8187 

0.7971 

5.7225 

5.4298 

9.930 

9.666 

2 

30 

0.44 

9495 

0.9958 

1.0023 

5.6546 

5.6489 

12.088 

12.168 

3 

60 

0.74 

28904 

0.7200 

0.7449 

5.6331 

5.7783 

8.7262 

9.0297 

■ 

120 

2.10 

86219 

0.6989 

0.7050 

5.3333 

5.4306 

8.4689 

8.5433 

5 

121 

2.13 

87422 

0.3857 

0.3944 

4.6381 

4.7817 

7.1029 

2.3049 

6 

240 

7.61 

157609 

0.4549 

0.4068 

6.0918 

5.9687 

7.8571 

2.5487 


Results and Discussion 

The flow computations were performed on a NASA MSFC Silicon Graphics Power Challenge 
cluster using 8 processors. The computational time for a typical calculation is 4.5x1 O' 5 CPU 


s/grid/step/processor. Usually, 2000-6000 iterations are required to achieve an approximate 
convergence. The radiation computations were performed on a NASA AMES Silicon Graphics 










































































































Origin 2000 cluster using 16 processors. The computational time for a typical mn is 5.0xl0' 5 CPU 
s/grid/step/direction/processor. In general, 20-40 iterations are used to achieve an approximately 
converged solution. 

Figures 7 shows the computed convective heat flux contours for the 30 s and 40 s trajectory 
points on the PPO at launch + 30 s trajectory. Different color scales are used for the plug base, the 
engine ramp, and the rest of the body surface such that the base-heating physics can be revealed. 
Unique three-dimensional base heat-flux contours are formed for each trajectory point due to the 
interaction of the reverse jet and the lateral wall jet with the pillowed plug base and the base-bleed. 
Contrasting to the approximately symmetric heat flux contours at 30 s (before PPO), the 
asymmetric heating pattern after PPO can be seen clearly at the 40 s heat flux contours, where the 
convective heat fluxes are much higher on the plug-base pillow and flex seals of the off-engine side 
(the right-hand side) than those of the on-engine side (the left-hand side), although the power level 
at 40 s is about half as that at 30 s. It is postulated that this asymmetric heating is caused by the 
reduced and asymmetric base-bleed and possibly a longer residence time for afterburning (due to 
reduced power level and reduced free stream Mach number). The convective heating decreases 
with time and becomes negligible at 200 s and 280 s due to a weak reverse jet weakened by the 
reduced power level and the high altitude plume attenuation (not shown). As a result of the reduced 
power level, there is negligible lateral plume jet impingement heating on the outer- and offset outer- 
bases at high altitudes, when compared with that computed at the full-power level. 3 A similar 
phenomenon occurs at the 60 s (before PPO) and 61 s (after PPO) for the PPO at launch + 60 s 
trajectory. For PPO at launch + 120 s trajectory, it appears that the PPO imposes negligible effect 
on base-heating due to low power level and high altitude. 

Figures 8 shows the computed radiative heat flux contours for the 30 s and 40 s trajectory points 



on the PPO at launch + 30 s trajectory. The radiative heating is transported through “view factors”, 
whereas the convective heating is transported through direct contact of the propulsive flow with the 
solid surface. Hence, the computed radiative heat flux contours look very different from those of 
the convection. For example, the upper and lower portions of the inner and outer bases show effect 
of radiative heating but not the middle section, since the nozzle plug blocks the radiation from the 
upper and lower plumes. In general, it is found that the computed radiative heat flux after PPO is 
lower than that before PPO. That is, the effect of reduced power level on radiative heat flux is 
bigger than that of the reduced base-bleed. 

Figures 9 shows computed two-dimensional centerline convective heat fluxes versus wetted 
distance in the horizontal direction on the various base surfaces of interest for the three PPO 
scenarios. For clarity, only those at selected trajectory points (before and after PPO) are shown, 
while the centerline heat flux at 0 s is shown as a baseline for comparison. It can be seen from Fig. 
9, for the PPO at launch + 30 s trajectory, that although the engine power level (PL) is reduced to 
about 50% due to the PPO event, the peak convective heating on the off-engine side plug-base at 40 
s is much higher (and closer to the middle) than that at 30 s and 100% PL, caused by an uneven 
split of the reduced base-bleed that sends only 2.919 lb/s to the off-engine side and 8.995 lb/s to the 
on-engine side. In fact, the convective heating is higher than that at sea level (0 s), which in turn is 
higher than that at 30 s. At a later time (not shown), e.g., at 100 s, although there is still asymmetric 
heating on the plug-base, the absolute magnitude of the heating is less than that prior to the PPO, 
due to higher plume attenuation. At 200 s and 280 s, the plume expansion is too strong to 
numerically stabilize the special base-bleed boundary condition and a fixed inlet has to be used. It 
can be safely assumed that the convective heating at those altitudes is negligible based on the 
results from other base surfaces. At 61 s (1 s after PPO) of the PPO at launch + 60 s trajectory. 



there is asymmetric heating and the heating levels are higher than those of the 60 s (before PPO) on 
both on- and off-engine side. The heating level at 60 s is already low due to already reduced power 
level. Therefore, the heating level at 61 s is lower than that at 0 s and 30 s. Hence, it is determined 
that PPO has negligible impact on plug-base at PPO + 60 s trajectory. Similarly, it is concluded 
that PPO has negligible impact on plug-base heating to the PPO + 120 s trajectory as well. 

Figure 10 shows the effect of PPO on the lower flex seal horizontal centerline convective heat 
fluxes for the three PPO scenarios. The heating profiles for the upper and lower seals are slightly 
different due to the asymmetry of the upper and lower geometries of the X-33 vehicle. The heating 
levels for the two seals are very similar, so only those on the lower seals are discussed. It can be 
seen that the curves in PPO at 30 s trajectory have distinctive peaks at the center, caused by the 
strong wall jet formed by the reverse jet impingement between the pillows. This strong wall jet has 
no place to go but upwards (toward the upper flex seal) and downwards (toward the lower flex 
seal), where it is then swiped away by the propulsive flow to form the double triangle recirculation 
pattern. 4 The two “wings” on both sides of the central heating curve are caused by the wall jet 
formed by the reverse jet on the pillows, which swirls down the pillow and flows through the seal 
region. The central heating level of the 30 s curve is higher then that of the 0 s (due to higher 
engine power level), but its “wing” heating is flattened out and lower than that of the 0 s (due to 
higher altitude and more base-bleed). At 40 s, although the central heating level is less than that of 
the 30 s (due to reduced power level during PPO), its “wing’ heating level is higher and closer to 
the center (due to reduced base-bleed); and it is higher at the “off” engine side (due to the 75/25 
split). The curve at 61 s (after PPO) is higher than that at 60 s (before PPO) in PPO at L + 60 s 
trajectory. The curves in the PPO at 120 s trajectory are both low and the effect of PPO is 
negligible. The heating for the PPO at launch + 60 s and PPO at launch + 120 s trajectories is less 


than that at launch + 30 s since base convective heating is less significant at higher altitudes. 

Figures 11 and 12 show radiative heat fluxes similar to those shown in Fig’s 9-10 for convective 
heat fluxes for the three PPO trajectories. In addition to the view factor, the radiative heating levels 
are also functions of power level, altitude, base-bleed, and ffee-stream Mach number. Of all the 
cases studied, the radiative heating after PPO is mostly lower than that seen before PPO on all base 
surfaces. This is particularly true for the radiative heating on the flex seals. As altitude or PPO 
time increases, the radiative heating decreases. Note that the level of radiative heating is quite 
different from that of convective heating on various components. For example, the peak level of 
radiative heating (Fig. 11) on the base horizontal centerline for PPO at launch + 30 s trajectory is 
about twice as that of convective heating (Fig. 9), whereas the peak heating levels on flex seals for 
both (Fig’s 10 and 12) are about the same. Again, the PPO has no effect on the PPO at L + 120 s 
trajectory since the heating level before PPO is already negligible. 

Conclusions 

A computational heat transfer analysis was performed to study the linear aerospike plume 
induced base-heating environment during PPO in ascent flight. Eighteen trajectory points selected 
from three PPO ascent scenarios were computed. The results are presented as three-dimensional 
convective and radiative heat flux contours and two-dimensional convective and radiative heat flux 
profiles versus wetted distance on component surfaces of interest. The base-heating environment 
during PPO is enveloped and the dual-engine base-heating physics during PPO is illustrated. In 
general, the locations that endure the most severe heating before and after PPO are the center of the 
flex seal and the center of the plug-base. PPO alters convective base-heating most when it occurs 
in the earliest scenario (PPO at launch + 30 s trajectory). The source of the impact of PPO on base- 


heating comes from the asymmetric and reduced base-bleed. 
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Fig. 1 Layout of an X-33 surface computational grid with component definitions. 
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Fig. 5 A comparison of radiative heat fluxes for radiometer R105. 
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Fig. 6 The thrust, altitude, Mach number profiles and computational points. 





Fig. 7 Convective heat flux contours for power-pack out at launch + 30 s trajectory. Color scales: 
0 - 400 for plug base, 0 - 9000 for engine ramp, and 0 - 60 for the rest of the body surface. 



Fig. 8 Radiative heat flux contours for power-pack out at launch + 30 s trajectory. Color scales: 0 
- 600 for plug base, 0 - 600 for engine ramp, and 0 - 170 for the rest of the body surface. 
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Fig. 9 Effect of power-pack out on base horizontal centerline convective heat fluxes for three 
power-pack out trajectories. 
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Fig. 1 1 Effect of power-pack out on base horizontal centerline radiative heat fluxes for three 
power-pack out trajectories. 
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